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The spectral distortion in the Cosmic Microwave Background Radiation (CMBR) caused from 
the inverse Compton scattering of low energy CMB photons by thermally distributed hot electrons 
in massive clusters of galaxies - better known as the thermal Sunyaev-Zel'dovich effect is usually 
estimated by the analytical or semi-analytical solutions of Kompaneets equation. These solutions 
are based on the assumption that the scattering is isotropic, i.e., the intensity of the CMBR remains 
angle-independent after scattering. In this letter, I present the solution of the full radiative transfer 
equations by incorporating anisotropic inverse Compton scattering of CMBR photons by relativistic 
electrons and show that the thermal Sunyaev-Zel'dovich effect is strongly angle-dependent. The 
spectral distortion in the CMBR alters significantly in different angular directions. The amount of 
distortion increases several times when the angle between the axis of symmetry and the photon ray 
path is increased. It is pointed out that for anisotropic scattering the Compton y parameter is angle 
dependent. Hence it is important to consider the direction of the emergent radiation with respect 
to the axis of symmetry in the electron rest frame while estimating any cosmological quantity. 
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Inverse Compton scattering of the Cosmic Microwave Background Radiation (CMBR) by hot intra-cluster gas - 
better known as the thermal Sunyaev-Zel'dovich effect (TSZE) [1,2] results in a systematic transfer of photons from 
the Raylcigh- Jeans to the Wien side of the spectrum causing a distortion in the Planckian nature of the spectrum. 
The measurements of the effect yield directly the properties of the hot intra-cluster gas, the total dynamical mass of 
the cluster as well as the indirect information on the cosmological evolution of the clusters. The effect is also used as 
an important tool to determine the Hubble constant Ho and the density parameter fio of the Universe [3-9]. Recent 
interferometric imaging of the TSZE [10,11] has been used to estimate the mass of various galaxy clusters which could 
constrain the cosmological parameters of structure formation models. Measurements of TSZE and the kinetic Sunyaev 
Zel'dovich effect (KSZE) from the Sunyaev-Zel'dovich Infrared Experiment and OVRO/BIMA determine the central 
Compton y parameter and constrain the radial peculiar velocity of the clusters. In near future sensitive observations 
of the effect with ground based and balloon-borne telescopes, equipped with bolometric multi-frequency arrays, are 
expected to yield high quality measurements. However, in modeling the observational data, the analytical solution 
of Kompaneets equation [12] with relativistic corrections are used in general. The Kompaneets equation describes 
the isotropic intensity of the radiation field after scattering. This means the expression used so far does not describe 
the anisotropy or the angle-dependence of the radiation after scattering. Hence the angle dependence of the spectral 
distortion in CMBR is completely ignored in any theoretical discussion or interpretation of the observed data in spite 
of the fact that scattering of an initially isotropic radiation should become anisotropic. In this letter I show that the 
thermal Sunyaev-Zel'dovich effect is strongly angle dependent. 

Chandrasekhar [13] has provided the radiative transfer equations for the anisotropic Compton scattering by assum- 
ing that the electron energy is much less than the photon energy. On the other hand, for the scattering of CMBR, 
the photon energy is much less than the electron energy. In fact, in the TSZE, the electrons are considered to have 
relativistic motion described by relativistic Maxwellian distribution. The relevant radiative transfer equations that 
describe the anisotropic Compton scattering of low energy photons can be written as [14-16] 
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Here ujq is the albedo for single scattering, \x = cos£? where 9 is the angle between the axis of symmetry (z axis) 
in the rest frame of electron gas and the ray path, t{z) is the total optical depth along the axis of symmetry given 
by rfr(z) = J <jTn e {z)dz, ot and n e being the Thomson scattering cross section and the electron number density 
respectively, k, c, m e , v and T e are Boltzmann constant, velocity of light, electron rest mass, frequency of the photon 
and the temperature of the electron respectively, ujq = 1 for a purely scattering medium. The phase functions Po, Pi 
and Pi are given as : 



P (m,m') = g [3 



H (1 



3m 2 )] 



(2) 



1 



and 



Pi(/i, //) = i [1 - 3/i' 2 - 3/i 2 (l - 3/x' 2 ) + 2/i 3 //(3 - 5 M ' 2 ) + 2/x//(3// 2 - 1)] , (3) 



P 2 (,x, M ') = | [3 - /i 2 - // 2 + nn'^/ip' - 5 + 3 M 2 + 3// 2 - 5 M 2 // 2 )] • (4) 



For an isotropic radiation field ^ J_ 1 I(fj,' ,v)dfj.' = I{v), \ f_- i I(fi' ,v)[i! 2 d[i! = and J_ x 7(ju' ,v)fi' 'd/j' = 

f \ I(fi, is)n' 3 dn' = 0.0 Therefore, for isotropic scattering with s as the ray path that does not change after scattering, 
equation (1) reduces to 

57(i/) _ 2kT e ( t z &I(y) _ 2v dlWY\ (5) 



dr(s) m e c 2 \ dv 2 dv 

which is well known as the Kompaneets equation [12] for low frequency. Note that for plane parallel medium, the 
operator d/ds becomes fxd/dz [13] which incorporates the change in the ray path with respect to the z axis after 
scattering. 

The analytical solution of equation (5) provides the thermal component of the distortion AI and is written as 
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where x — kT ^ IBR , h — 

2(kT CM BB.) 3 /(hc) 2 . The term y=^J <jTn e ds is usually referred to as the Compton y 
parameter and St is a relativistic correction to the thermal effect [17] significant if kT e > lOKeV. 

The above expression is usually used in order to estimate the TSZE. Clearly, the angle dependence of the emergent 
intensity and hence the distortion is neglected as the radiation field after scattering is assumed to remain isotropic. 

Now equation (1) is a coupled integro partial differential equation and so it cannot be solved analytically. I solve it 
numerically by discretization method. In this method the medium is divided into several shells and the integration is 
performed over two dimensional grids of angular and radial points. For the angle integration I have adopted an eight 
point Gauss-Legendre roots and weights. I have taken eighty frequency points with equal spacing. For the initial 
condition, I have provided equal amount of intensity corresponding to Tqmbr = 2.728K at r = along all directions 
i.e., 7(/i, v,t = 0) = ^(e 1 — The code is thoroughly tested for stability and flux conservation. The numerical 

results coincides with the analytical solution given by equation (6) when the radiation is made isotropic. On the other 
hand, if P\ = P^ = 0, the results are well matched with that presented in [13] for Thomson scattering. 

Usually the isothermal (3 models [18] is considered for the density distribution of the clusters. The spherical 
isothermal model density is described by 

/ r 2 \ -3/3/2 

n e (r) = no I 1 + -J (7) 



where the core radius r c and j3 are shape parameters, uq is the central electron density. However, for the present 
purpose it is sufficient to consider an isothermal, homogeneous and plane parallel medium with a constant optical 
depth. In the present work I have taken a constant value of n e = 10~ 3 and the size of the cluster is 4 Mpc. The 
electron temperature T e is taken to be 1 KeV and 10 KeV. The results are presented graphically in Fig. 1 and in Fig. 
2. 

The spectral distortion for the isotropic case is characterized by three distinct frequencies : the crossover frequency 
xq = 3.83 where the TSZE vanishes; x m i n = 2.26 which gives the minimum decrement of the CMB intensity and 
Xmax = 6.51 which gives the maximum distortion due to this effect. The value of xo is however pushed to higher 
values of x with the increase in T e for the relativistic case. 

First of all there is no change in the values of xq, x m %n and x max in any direction even if the radiation field is 
anisotropic. In order to show this, the result for the isotropic scattering is also presented in Fig. 1 and in Fig. 2. 
I have taken the ray path for the isotropic case along the axis of symmetry. The values of the three characteristic 
frequency points remain unchanged because of the fact that the frequency dependence of the distortion is independent 
of the angle dependence of the distortion as can be seen from equation (1). 



l In Ref. [15] I = \l should be 7 = \l . 
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FIG. 1. Anisotropic Sunyaev-Zel'dovich effect with kT e — 1 KeV. From top to bottom the solid lines represent the spectral 
distortion due to anisotropic scattering along the angular directions = O.f, (2) fi — 0.24, (3) fi = 0.4, (4) fi = 0.6, 

(5)/U = 0.76, (6) /i = 0.9, (7) fi = 0.98. The dashed line represent the distortion due to isotropic scattering with the ray path 
being along the axis of symmetry. 
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FIG. 2. Same as Fig.f but with kT e = 10 KeV. 
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It is worth mentioning at this point that the photon changes its direction each time it gets scattered. With respect 
to the z axis, the new ray path becomes s — z/fi. As a consequence, the optical depth through the new ray path 
becomes dr(z)/ ' \i where dr(z) corresponds to the optical depth along the z axis. Therefore, for anisotropic scattering 
it is incorrect to take dl/dr(s) instead of ndl/dr(z) in the left hand side of equation (1) because dl/dr(s) describes 
the change in the intensity along the same ray path s before and after scattering. Flux conservation confirms this fact. 
The Compton y parameter, for a single anisotropic scattering in a plane-parallel medium, therefore, can be defined as 



where fi is now the angle between the incident direction to the emergent angle of the radiation. In reality, however, 
the photon would suffer multiple scattering depending on the density of the medium and hence would change its 
direction several times before it emerges out. If the medium is spherically symmetric then the situation becomes more 
complicated. In plane parallel geometry the ray makes a constant angle 9 with the normal while in the spherically 
symmetric geometry the angle made by the ray direction and the radius vector changes constantly. 

Now, the most interesting and important message conveyed by Fig.l and Fig. 2 is that the spectral distortion is 
strongly angle dependent. The degree of distortion in the CMBR due to TSZE changes drastically with the change in 
the angular direction of the photons. This clearly demonstrates the fact that the TSZE produces anisotropic distortion 
in the CMBR spectrum. 

The degree of distortion increases as the angle between the emergent intensity and the axis of symmetry increases. 
Therefore the minimum decrement and maximum distortion of CMBR occur at x = 2.26 and x = 6.51 respectively 
for the radiation emerging almost perpendicular to the axis of symmetry. For a single scattering this can be explained 
easily. With the increase in 9 and hence with the decrease in n the Compton y parameter increases. As a result the 
distortion increases when 9 increases. For a plane parallel medium the Compton y parameter becomes infinite when 
9 = tt/2 and hence the distortion becomes infinite at 9 — tt/2. It should be mentioned here that the results presented 
in Fig. 1 and in Fig. 2 incorporate multiple scattering of photons incident and emergent at any angular direction 
before and after scattering. The radiation traversing at the opposite direction, i.e., the backscattered radiation is also 
taken into care while solving equation (1) numerically. The qualitative nature of the result will not alter if we consider 
spherical symmetry except the fact that the ray would peak with the radius vector towards the outer boundary of the 
sphere. The results with spherical symmetry will be published in a forthcoming paper. 

Fig. 1 and Fig. 2 show that the spectral distortion due to isotropic scattering matches with that due to anisotropic 
scattering when fx = 0.76. Therfore the distortion due to anisotropic scattering would be less as compared to that due 
to isotropic scattering if fi < 0.76. Now, the spectral distortion along large angular direction should be irrelevant from 
the observational point of view. But the important point to be noted is that- if the effect is measured along the line 
of sight ( which can be assumed as the axis of symmetry) then for the same value of the Compton y parameter, the 
distortion due to anisotropic scattering would be much less than that calculated under isotropic assumption. In other 
words, the spectral distortion along the line of sight would result into underestimation of the Compton y parameter 
if isotrpic scattering is assumed. 

Therefore it is extremely important to know the direction of the observed intensity with respect to the axis of 
symmetry in the rest frame of the electron gas while estimating the density of the clusters. The isotropic assumption 
would either underestimate or overestimate it depending on the angular direction of the emergent intensity of CMBR. 

In conclusion I would like to emphasize that the anisotropy in the CMBR, induced by Compton scattering from hot 
electrons in the intra-cluster, not only yields anisotropic distortion in the Planckian spectrum but also would results 
into polarization in the CMBR. 

I am thankful to P. Bhattacharjee for bringing my attention on this problem and for discussions. 
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